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Abstract 

Various new physics models predict a light CP-odd Higgs boson (labeled as a) and open up new 
decay modes for Z-boson, such as Z — ffa, Z — ay and Z — aaa, which could be explored at 
the GigaZ option of the ILC. In this work we investigate these rare decays in several new physics 
models, namely the type-II two Higgs doublet model (type-II 2HDM), the lepton-specific two 
Higgs doublet model (L2HDM), the nearly minimal supersymetric standard model (nMSSM) and 
the next-to-minimal supersymmetric standard model (NMSSM). We find that in the parameter 
space allowed by current experiments, the branching ratios can reach 1074 for Z > ffa (f = b,7), 
10-9 for Z + ay and 107? for Z — aaa, which implies that the decays Z — f fa and Z — aaa 
may be accessible at the GigaZ option. Moreover, since different models predict different patterns 
of the branching ratios, the measurement of these rare decays at the GigaZ may be utilized to 


distinguish the models. 


PACS numbers: 13.38.Dg,12.60.Fr,14.80.Da 


I. INTRODUCTION 


The LEP experiments at the resonance of Z-boson have tested the Standard Model (SM) 
at quantum level, measuring the Z-decay into fermion pairs with an accuracy of one part in 
ten thousands. The good agreement of the LEP data with the SM predictions have severely 
constrained the behavior of new physics at the Z-pole. Taking these achievements into 
account one can imagine that the physics of Z-boson will again play the central role in the 
frontier of particle physics if the next generation Z factory comes true with the generated 
Z events several orders of magnitude higher than that of the LEP. This factory can be 
realized in the GigaZ option of the International Linear Collider mon The ILC is a 
proposed electron-positron collider with tunable energy ranging from 400GeV to 500GeV 
and polarized beams in its first phase, and the GigaZ option corresponds to its operation 
on top of the resonance of Z boson by adding a bypass to its main beam line. Given the 
high luminosity, L = 7 x 10cm ?s^!, and the cross section at the resonance of Z boson, 
cz œ 30nb, about 2 x 10? Z events can be generated in an operational year of 107s of 
GigaZ, which implies that the expected sensitivity to the branching ratio of Z-decay can 
be improved from 10^? at the LEP to 1078 at the GigaZ il. In light of this, the Z-boson 
properties, especially its exotic or rare decays which are widely believed to be sensitive to 
new physics, should be investigated comprehensively to evaluate their potential in probing 
new physics. 


studied to explore the flavor texture in new physics |2], and it was found that, although these 


Among the rare Z-decays, the flavor ate los processes were most extensively 
processes are severely suppressed in the SM, their branching ratios in new physics models 
can be greatly enhanced to 1078 for lepton flavor violation decays and 10-9 for quark 
flavor violation decays |4]. Besides the FC processes, the Z-decay into light Higgs boson(s) 
is another type of rare process that was widely studied, e.g. the decay Z > f fa (f = 6,7) 
with the particle a denoting a light Higgs boson was studied in n the decay Z — ay 
was studied in the two H i” doublet model (2HDM)|6] and the minimal supersymmetric 
standard model (MSS , and the decay Z — aaa was studied in a model independent 
way B. in veal ha sl and also in MSSMIiI. These studies indicate that, in contrast 
with the kinematic forbidden of these decays in the SM, the rates of these decays can be 


as large as 107? in new physics models, which lie within the expected sensitivity of the 


GigaZ. In this work, we extend the previous studies of these decays to some new models 
and investigate these decays altogether. We are motivated by some recent studies on the 
singlet extension of the MSSM, such as the next-to-minimal supersymmetric standard model 
(NMSSM) 13 and the nearly minimal supersymmetric standard model (nMSSM) ' 
where a light CP-odd Higgs boson a with singlet-dominant component may naturally arise 
from the spontaneous breaking of some approximate global symmetry like Up(1) or Peccei- 
Quuin symmetry 17]. These non-minimal supersymmetric models can not only avoid 
the u-problem, but also alleviate the little hierarchy by having such a light Higgs boson 
a |18]. We are also motivated by that, with the latest experiments, the properties of the 
light Higgs boson are more stringently constrained than before. So it is worth updating the 
previous studies. 

5o far there is no model-independent lower bound on the lightest Higgs boson mass. In 
the SM, it must be heavier than 114 GeV, obtained from the null observation of the Higgs 
boson at LEP experiments. However, due to the more complex structure of the Higgs sector 
in the extensions of the SM, this lower bound can be significantly relaxed according to recent 
studies, e.g., for the CP-odd Higgs boson a we have m,21 GeV in the nMSSM Maz 0.21 
GeV in the NMSSM |16], and m,27 GeV in the lepton-specific 2HDM (L2HDM) "i With 
such a light CP-odd Higgs boson, the Z-decay into one or more a is open up. Noting that the 
decay Z — aa is forbidden due to Bose symmetry, we in this work study the rare Z-decays 
Z — ffa (f —b,7T), Z ^ ay and Z — aaa in a comparative way for four models, namely 
the Type-II zu the L2HDM hd, the nMSSM and the NMSSM. In our study, 
we examine carefully the constraints on the light a from many latest experimental results. 

This work is organized as follows. In Sec. II we briefly describe the four new physics 
models. In Sec. III we present the calculations of the rare Z-decays. In Sec. IV we list the 
constraints on the four new physics models. In Sec. V we show the numerical results for the 
branching ratios of the rare Z-decays in various models. Finally, the conclusion is given in 


Sec. VI. 


I. THE NEW PHYSICS MODELS 


As the most economical way, the SM utilizes one Higgs doublet to break the electroweak 
symmetry. As a result, the SM predicts only one physical Higgs boson with its properties 


totally determined by two free parameters. In new physics models, the Higgs sector is usually 


extended by adding Higgs doublets and/or singlets, and consequently, more physical Higgs 
bosons are predicted along with more free parameters involved in. 
The general 2HDM contains two SU(2); doublet Higgs fields ¢, and and with the 


assumption of CP-conserving, its scalar potential can be parameterized as|21 


V = miélós + molda — (mólos + H.c) + 2 (01) + Zhe)? 
A 
«(61i (662) + A (612) (0501) + 77 (0102)? + He], (1) 


where A; (i = 1,--- ,5) are free dimensionless parameters, and m; (i = 1,2,3) are the 
parameters with mass dimension. After the electroweak symmetry breaking, the spectrum 


of this Higgs sector includes three massless Goldstone modes, which become the longitudinal 


modes of W* and Z bosons, and five massive physical states: two CP-even Higgs bosons hı 


and A5, one neutral CP-odd Higgs particle a and a pair of charged Higgs bosons H^. Noting 
the constraint v? J- v? = (246 GeV)? with v; and v9 denoting the vacuum expectation values 


(vev) of $1 and $» respectively, we choose 


Mh, Mh, Ma, mgr, tanp, sina, As (2) 


as the input parameters with tan 3 = v2/v,, and a being the mixing angle that diagonalizes 
the mass matrix of the CP-even Higgs fields. 

The difference between the Type-II 2HDM and the L2HDM comes from the Yukawa 
coupling of the Higgs bosons to quark/lepton. In the Type-II 2HDM, one Higgs doublet $» 
generates the masses of up-type quarks and the other doublet o, generates the masses of 
down-type quarks and charged leptons; while in the L2HDM one Higgs doublet $4 couples 
only to leptons and the other doublet $5 couples "ubl So the Yukawa interactions 


of a to fermions in these two models are given by |20, [21 


igm 


1gMu, igm 


Dea = 27, 29 bui uat > tan pdi dia + ,— tan Bey eia (3) 
igm,, j 

Dye = 197773 cot būiy uia — L s cot Bdyy?dya + —— tom tan Bē eia (4) 
2mw 2mw 2m 


with i denoting generation index. Obviously, in the Type-II 2HDM the bba coupling and 
the Tra coupling can be simultaneously enhanced by tan 3, while in the L2HDM only the 


TTG coupling is enhanced by tan f. 


The structures of the nMSSM and the NMSSM are described by their superpotentials 


and corresponding soft-breaking terms, which are given by |24 


Wamssm = Wmussm + AH, - Ha + £ M26, (5) 
Wyussm = Wussm + AH, - HaS + Engs, (6) 
VIROM = m?IH.? + mH? + MIS] + (AAS H, Hat Es M?S +h.c.), (7) 
VSM = m2|H,|? + MAHA + mzsp + (AASH, -Ha + id +h.c.), (8) 


where Wmssm is the superpotential of the MSSM without the u term, Haa and $ are Higgs 
doublet and singlet superfields with H,,; and S being their scalar component respectively, 
Mu, Ma, Ms, Ax, A, and ég M} are soft breaking parameters, and À and « are coefficients 
of the Higgs self interactions. 

With the superpotentials and the soft-breaking terms, one can get the Higgs potentials of 
the nMSSM and the NMSSM respectively. Like the 2HDM, the Higgs bosons with same CP 
property will mix and the mass eigenstates are obtained by diagonalizing the corresponding 


mass matrices: 


hy Pu a Pu 
H+ H+ 

hə U Qa ’ A =U Yd ’ =U ’ (9) 
Gr Hj 

ha o G? E 


where the fields on the right hands of the equations are component fields of H,, Hq and S 
defined by 


tat dated Ht 1 
H4 = : E H, = s D S = —= (s +0 + i$), (10) 
Hr — /2 


hı, h5, h3 and a, A are respectively the CP-even and CP-odd neutral Higgs bosons, G? and 
Gt are Goldstone bosons eaten by Z and W*, and H* is the charged Higgs boson. So both 
the nMSSM and NMSSM predict three CP-even Higgs bosons, two CP-odd Higgs bosons 
and one pair of charged Higgs bosons. In general, the lighter CP-odd Higgs a in these model 
is the mixture of the singlet field € and the doublet field combination, cos by, 4- sin bya, i.e. 


a = cos 04€ + sin 04(cos By, + sin bya), (11) 


Tu tan BsinO4. So for 


and its couplings to down-type quarks are then proportional to "m 


singlet dominated a, sin 04 is small and the couplings are suppressed. As a comparison, the 
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interactions of a with the squarks are given byld 


Lağ = pam (Av cot B cos 04 + A, cot B sin4 — usin 04) (Ahù — uj üg)a 
igmaq 


= (Av tan 8cosĝa + Agtan B sin 04 — usin 04) (dhdi — d;dg)a, (12) 
w 


ie. the interaction does not vanish when sin@4 approaches zero. Just like the 2HDM 
where we use the vevs of the Higgs fields as fundamental parameters, we choose A, &, tan f, 
Heff = A(S), A, and mà = s (Aa + *£) as input parameters for the NMSSM[24 and 
r, tan B, uer; = A(S), Ay, Ms and m3 = suas (HA, + àr M?) as input parameters for the 
nMSSM[I9|. 


About the nMSSM and the NMSSM, three points should be noted. The first is for the 


two models, there is no explicit jj—term, and the effective u parameter (uer) is generated 
when the scalar component of $ develops a vev. The second is, the nMSSM is actually same 
as the NMSSM with « = olg. because the tadpole terms £e M2$ and its soft breaking 
term £sM?S in the nMSSM do not induce any interactions, except for the tree-level Higgs 
boson masses and the minimization conditions. And the last is despite of the similarities, 
the nMSSM has its own peculiarity, which comes from its neutralino sector. In the basis 


(=i), 3,0 We Vs), its neutralino mass matrix is given by [14 


Mi 0 Mzswsp —MzSw Cb 0 
0 Mo —mzcwsy MzZCwCb 0 
mzgSwsSy —MzcCwSp 0 -u — AUC, (13) 
—MzswCh mzCwCb —u 0 — AUS, 
0 0 —ÀAUCh — AUS, 0 


where M; and M» are U(1) and SU(2) gaugino masses respectively, sw = sinÜw, cw = 
cos Ow, sy, = sin and cy = cos. After diagonalizing this matrix one can get the mass 
eigenstate of the lightest neutralino X? with mass taking the following form " 


2uA? (v2 + v2) tan 8 


ot —Ó——— 14 
2u? + M(v2 + v2) tan? B+1 9) 


m 


0 C 
1 


This expression implies that X? must be lighter than about 60 GeV for > 100GeV (from 
lower bound on chargnio mass) and A < 0.7 (perturbativity bound). Like the other super- 
symmetric models, X? as the lightest sparticle acts as the dark matter in the universe, but 


due to its singlino-dominated nature, it is difficult to annihilate sufficiently to get the correct 
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density in the current universe. So the relic density of X? plays a crucial way in selecting 
the model parameters. For example, as shown in ld. for X? > 37GeV, there is no way to 
get the correct relic density, and for the other cases, X? mainly annihilates by exchanging 
Z boson for 30GeV < myo < 37GeV , or by exchanging a light CP-odd Higgs boson a with 
mass satisfying the relation Ma ~ 2myo for mgo < 25GeV. For the annihilation, tan and 
u are required to be less than 10 and 500GeV respectively because through Eq. (14) a large 
tan B or u will suppress m zo to make the annihilation more difficult. The properties of the 
lightest CP-odd Higgs boson a, such as its mass and couplings, are also limited tightly since 


a plays an important role in X? annihilation. ul e of the nMSSM is also 


rather special, and this was discussed in detail in [19]. 


IIl. CALCULATIONS 


T/b T/b -a 
Z 
4 
xv ^ d 
5 
NSS 
Z a Z "E Z Re T/b 
h 
7/b 7/b 7/b 


(a) (b) (c) 
FIG. 1: Feynman diagrams contributing to the decay Z > f fa (f = b,7) in new physics models. 


h denotes all possible intermediate CP-even Higgs bosons in the corresponding model. 


FIG. 2: Feynman diagrams contributing to the decay Z — ay at one-loop level in new physics 


models. Note the chargino loop diagrams only exist in the nMSSM and NMSSM. 


In the Type-II 2HDM, L2HDM, nMSSM and NMSSM, the rare Z-decays Z — ffa 
(f —b,T), Z > ay and Z — aaa may proceed by the Feynman diagrams shown in Fig [i] 
Fig[2] and Fig[3] respectively. For these diagrams, the intermediate state h represents all 


T 


avavavavas q Yee NK, " " e Yen i em 
Z > 2a Z h e Z h A 
h^ p x ih xc jd 
x ay hrs 
N N N. 
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Ra MEN k k 
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FIG. 3: Feynman diagrams contributing to the decay Z — aaa in new physics models. Only the 
correction from the Higgs-mediated loops are considered since the other corrections can be safely 


neglected (see the arguments in the text). 


possible CP-even Higgs bosons in the corresponding model, i.e. ^; and hg in Type-II 2HDM 
and L2HDM and hy, h3 and h3 in nMSSM and NMSSM. In order to take into account the 
possible resonance effects of h in Fig[I{c) for Z — f fa and FigB] (a) for Z — aaa, we have 
calculated all the decay modes of h and properly included the width effect in its propagator. 
As to the decay Z — ay, two points should be noted. One is, unlike the decays Z > ffa 
and Z — aaa, this process proceeds only through loops mediated by quarks/leptons in the 
Type-II 2HDM and L2HDM, and additionally by sparticles in the nMSSM and NMSSM. 
5o in most cases its rate should be much smaller than the other two. The other is due to 
CP-invariance, loops mediated by squarks/sleptons give no contribution to the decay|7]. In 
actual calculation, this is reflected by the fact that the coupling coefficient of gpqza differs 
from that of d; ága by a minus sign (see Eq.(12)), and as a result, the squark-mediated 
contributions to Z — ay are completely canceled out. 

With regard to the rare decay Z — aaa, we have more explanations. In the lowest 
order, this decay proceeds by the diagram shown in Fig[3] (a), and hence one may think 
that, as a rough estimate, it is enough to only consider the contributions from Fig[3[a). 
However, we note that in some cases of the Type-II 2HDM and L2HDM, due to the can- 
celation of the contributions from different h in Fig] (a) and also due to the potentially 


largeness of haa couplings (i.e. larger than the electroweak scale v), the radiative correc- 
tion from the Higgs-mediated loops may dominate over the tree level contribution even 
when the tree level prediction of the rate, Bripee(Z — aaa), exceeds 10 9. On the other 
hand, we find the contribution from quark/lepton-mediated loops can be safely neglected if 
Briree(Z — aaa) > 107? in the Type-II 2HDM and the L2HDM. In the nMSSM and the 
NMSSM, besides the corrections from the Higgs- and quark/lepton-mediated loops, loops 
involving sparticles such as squarks, charginos and neutralinos can also contribute to the 
decay. We numerically checked that the contributions from squarks and charginos can be 
safely neglected if Bripee(Z — aaa) > 1078. We also calculated part of potentially large 
neutralino correction (note that there are totally about 5* diagrams for such correction!) 
and found they can be neglected too. Since considering all the radiative corrections will 
make our numerical calculation rather slow, we only include the most important correction, 
namely that from Higgs-mediated loops, in presenting our results for the four models. 

One can intuitively understand the relative smallness of the sparticle contribution to 
Z — aaa as follows. First consider the squark contribution which is induced by the Zq*q 
interaction (d denotes the squark in chirality state) and the q*ga interaction through box 
diagrams. Because the Zq*q interaction conserves the chirality of the squarks while the 
d'qa interaction violates the chirality, to get non-zero contribution to Z — aaa from the 
squark loops, at least four chiral flippings are needed, with three of them provided by q*qa 
interaction and the rest provided by the left-right squark mixing. This means that, if one 
calculates the amplitude in the chirality basis with the mass insertion method, the amplitude 
is suppressed by the mixing factor Taya with m,Xq being the off diagonal element in squark 
mass matrix. Next consider the etis rettet contributions. Since for a light a, 
its doublet component, parameterized by sin 04 in Eq.(11), is usually small, the couplings 
of a with the sparticles will never be tremendously imr) So the chargino/neutralino 
contributions are not important too. In our calculation of the decays, we work in the mass 


eigenstates of sparticles instead of in the chirality basis. 


IV. CONSTRAINTS ON THE NEW PHYSICS MODELS 


For the Type-II 2HDM and the L2HDM, we consider the following constraints f 


(1) Theoretical constraints on A; from perturbativity, unitarity and requirements that the 


scalar potential is finit at large field values and contains no flat directions lal bd. 
which imply that 


eae uml. Aue, Aso Lo Ais As eA; 

3( + Ao) + / 90 — A2)? + 4(2A3 + A4)? < 167, 

M+ Ag E/O — Ag)? + ADAGP « 167, 

Mt A2 E VA — A3)? + AD? < 167, 

As +24 + 3]A5s| < 81, A3 A4 < 81, A3 |As| < 8m. (15) 


(2) The constraints from the LEP search for neutral Higgs bosons. We compute the 
signals from the Higgs-strahlung production e*e^ — Zh; (i = 1,2) with h; > 
2b, 27, 4b, 47, 2027 
h;a — 4b, 47, 2b27, 6b, 67 30 , and compare them with the corresponding LEP data 


and from the associated production e*e^ — h;a with 


which have been inputted into our code. We also consider the constraints from 
ete” — Zh, by looking for a peak of Mj, recoil mass distribution of Z-boson [31 


and the constraint of P(Z — h;a) < 5.8 MeV when m, + mp, < mz |32]. 


These constraints limit the quantities such as C2/, = [977,,/9775,,,] X Br(hi > bb) on 
the C2 f — ma, plane with the the subscript gzzn, denoting the coupling coefficient of 
the ZZh; interaction. They also impose a model-dependent lower bound on m;j,, e.g., 
70GeV for the Type-II 2HDM (from our scan results), 50GeV for the L2HDM]20|, and 
30GeV for the nMSSM m These bounds are significantly lower than that of the SM, 
i.e. 114GeV, partially because in new physics models, unconventional decay modes of 
hj such as h; — aa are open up. As to the nMSSM, another specific reason for allowing 
a significantly lighter CP-even Higgs boson is that the boson may be singlet-dominated 


in this model. 


With regard to the lightest CP-odd Higgs boson a, we checked that there is no lower 


bound on its mass so long as the Zh;a interaction is weak or h; is sufficiently heavy. 


(3) The constraints from the LEP search for a light Higgs boson via the Yukawa process 
ete > ffS with f = b, r and S denoting a scalar led These constraints can limit 


the f f.S coupling versus mg in new physics models. 
(4) The constraints from the CLEO-III limit on Br(Y(1S) > ay > r*7-*) and the latest 
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BaBar limits on Br(T(3S) — ay > Ttr *y,u* uy). These constraints will put very 
tight constraints on the abb coupling for m, < 9GeV. In our analysis, we use the 


results of Fig.8 in the second paper of to excluded the unfavored points. 


The constraints from Zr*7- couplings. Since the Higgs sector can give sizable higher 
order corrections to Z7*7~ couplings, we calculate them to one loop level and require 
the corrected Zr*T- couplings to lie within the 2e range of their fitted value. The 
SM predictions for the couplings at Z-pole are given by gp" = —0.03712 and g3” = 
—0.50127 , and the fitted values are given by —0.0366 + 0.00245 and —0.50204 + 
0.00064, respectively 4l We adopt the formula in to the 2HDM in our calculation. 


The constraints from 7 leptonic decay. We require the new physics correction to the 
branching ratio Br(T — ev-v,) to be in the range of —0.80% ~ 1.21% |36]. We use 
the formula in in our calculation. 


About the constraints (5) and (6), two points should be noted. One is all Higgs 
bosons are involved in the constraints by entering the self energy of 7 lepton, the 
ZTT vertex correction or the W Tv, vertex correction, and also the box diagrams for 
T — eDV, bs kd. Since the Yukawa couplings of the Higgs bosons to 7 lepton get 
enhanced by tan 6 and so do the corrections, tan must be upper bounded for given 
spectrum of the Higgs sector. Generally speaking, the lighter a is, the more tightly 
tan is limitedl2d Cd The other point is in the Type-II 2HDM, Aj, B-physics 
observables as well as Y decays discussed above can constraint the model in a tighter 
way than the constraints (5) and (6) since the Yukawa couplings of 7 lepton and b 
quark are simultaneously enhanced by tan 8. But for the L2HDM, because only the 
Yukawa couplings of 7 lepton get enhanced (see Eq), the constraints (5) and (6) are 


more important in limiting tan f. 


Indirect constraints from the precision electroweak observables such as py, sin? 6' pf and 
My, or their combinations e;(7 = 1, 2,3) . We require e; to be compatible with the 
LEP/SLD data at 9596 confidence levell34 . We also require new physics prediction 
of Ry = T(Z — bb)/T(Z — hadrons) is within the 2e range of its experimental 
value. The latest results for Ry, are Rọ” = 0.21629 + 0.00066 (measured value) and 
Re“ = 0.21578 (SM prediction) for m, = 173 GeV ag. In our code, we adopt the 
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formula for these observables presented in to the Type-II 2HDM and the L2HDM 


respectively. 


In calculating pe, sin? Oe rt and My, we note that these observables get dominant 
contributions from the self energies of the gauge bosons Z, W and y. Since there is no 
Zaa coupling or yaa coupling, a must be associated with the other Higgs bosons to 


contribute to the self energies. So by the UV convergence of these quantities, one can 


infer that, for the case of a light a and mp,, my > mz, these quantities depend on the 


m2 
spectrum of the Higgs sector in a way like In "uc at leading order, which implies that 


a light a can still survive the constraints md the precision electroweak observables 


given the splitting between mp, and mjr« is moderate|20J. 


The constraints from B physics observables such as the branching ratios for B > Xy, 
B, > wtp” and Bt — rtv,, and the mass differences AM; and AM,. We require 
their theoretical predications to agree with the corresponding experimental values at 


2g level. 


In the Type-II 2HDM and the L2HDM, only the charged Higgs boson contributes 


to these observables by loops, so one can expect that my+ versus tan is to be 


limited. Combined analysis of the limits in the Type-II 2HDM has been done by the 
CKMfitter Group, and the lower bound of my+ as a function of tan 8 was given in 
Fig.11 of [39]. This analysis indicates that my+ must be heavier than 316GeV at 95% 
C.L. regardless the value of tan 8. In this work, we use the results of Fig.11 in 


to exclude the unfavored points. As for the L2HDM, B physics actually can not put 
any constraints 40 because in this model the couplings of the charged Higgs boson to 
quarks are proportional to cot @ and in the case of large tan 6 which we are interested 
in, they are suppressed. In our analysis of the L2HDM, we impose the LEP bound on 
My+, Le. My+ > 92GeV [41]. 


The constraints from the muon anomalous magnetic moment a,. Now both the theo- 
retical prediction and the experimental measured value of a,, have reached a remarkable 
precision, but a significant deviation still exists: a7"? — a7 = (25.5 +8.0) x 107? [42]. 
In the 2HDM, a,, gets additional contributions from the one-loop diagrams induced 
by the Higgs bosons and also from the two-loop Barr-Zee diagrams mediated by a and 


h;|43|]. If the Higgs bosons are much heavier than u lepton mass, the contributions 
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(10) 


from the Barr-Zee diagrams are more important, and to efficiently alleviate the dis- 
crepancy of a,,, one needs a light a along with its enhanced couplings to u lepton and 
also to heavy fermions such as bottom quark and 7 lepton to push up the effects of the 
Barr-Zee diagramlt3l. The CP-even Higgs bosons are usually preferred to be heavy 


since their contributions to a, are negative. 


In the Type-II 2HDM, because tan f is tightly constrained by the process ete~ — bba 
at the LEPB3 and the T decay l17, the Barr-Zee diagram contribution is insufficient to 
enhance a, to 2e range around its measured M» So in our analysis, we require 
the Type-II 2HDM to explain a, at 30 level. While for the L2HDM, tan £ is less 
constrained compared with the Type-II 2HDM, and the Barr-Zee diagram involving 
the 7-loop is capable to push up greatly the theoretical prediction of d Therefore, 


we require the L2HDM to explain the discrepancy at 2o level. 


Unlike the other constraints discussed above, the a, constraint will put a two-sided 
bound on tan 6 since on the one hand, it needs a large tan 8 to enhance the Barr-Zee 
contribution, but on the other hand, too large tan 9 will result in an unacceptable 


large ap- 


Since this paper concentrates on a light a, the decay h; — aa is open up with a 
possible large decay width. We require the width of any Higgs boson to be smaller 
than its mass to avoid a too fat Higgs boson[ig. We checked that for the scenario 
characterized by mj, /m, > 3, the coefficient of h;aa interaction is usually larger than 


the electroweak scale v, and consequently a large decay width is resulted. 


For the nMSSM and NMSSM, the above constraints become more complicated because 


in these models, not only more Higgs bosons are involved in, but also sparticles enter the 


constraints. So it is not easy to understand some of the constraints intuitively. Take the 


process B — X;,y as an example. In the supersymmetric models, besides the charged Higgs 


contribution, chargino loops, gluino loops as well as neutralino loops also contribute to the 


process|45], and depending on the SUSY parameters, any of these contributions may become 


dominated over or be canceled by other contributions. As a result, although the charged 


Higgs affects the process in the same way as that in the Type-II 2HDM, charged Higgs as 
light as 130GeV is still allowed even for tan 8 50/46]. 
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Since among the constraints, a, is rather peculiar in that it needs new physics to explain 


the discrepancy between af"? and a?” , we discuss more about its dependence on SUSY 


SM 
H m 


parameters. In the nMSSM and the NMSSM, a, receives contributions from Higgs loops and 
neutralino/chargino loops. For the Higgs contribution, it is quite similar to that of the Type- 
II 2HDM except that more Higgs bosons are involved Eel For the neutralino/chargino 
contribution, in the light bino limit (i.e. Mı < Mo, u), it can be approximated by|48 


100GeV V? / u — A, cot 
ja, = 18tan f (m cs) 10-19 (16) 


mg 

for mg, c mg, = mj = 2M, with mj, being smuon mass. So combining the two contribu- 

tions together, one can learn that a light a along with large tan 6 and/or light smuon with 
moderate tan @ are favored to dilute the discrepancy. 

Because more parameters are involved in the constraints on the supersymmetric models, 


we consider following additional constraints to further limit their parameters: 


(a) Direct bounds on sparticle masses from the LEP1, the LEP2 and the Tevatron exper- 


iments |38]. 


(b) The LEP1 bound on invisible Z decay '(Z — ¥?X?) < 1.76 MeV; the LEP2 bound 
on neutralino production o(ete~ — X1X1) < 107? pb (i > 1) and o(ete™ — X2X2) < 
10-! pb (i, j > 1) 49). 


(c) Dark matter constraints from the WMAP relic density 0.0975 « QA? « 0.1213 Ed. 


Note that among the above constraints, the constraint (2) on Higgs sector and the con- 
straint (c) on neutralino sector are very important. This is because in the supersymmetric 
models, the SM-like Higgs is upper bounded by about 100GeV at tree level and by about 
140GeV at loop level, and that the relic density restricts the LSP annihilation cross section 
in a certain narrow range. 

In our analysis of the NMSSM, we calculate the constraints (3) and (5-7) by ourselves 
and utilize the code NMSSM Tools to implement the rest constraints. We also extend 
NMSSM Tools to the nMSSM to implement the constraints. For the extension, the most 
difficult thing we faced is how to adapt the code micrOMEGAs{5J to the nMSSM case. We 
solve this problem by noting the following facts: 


14 


e As we mentioned before, the nMSSM is actually same as the NMSSM with the trilinear 


singlet term setting to zero. So we can utilize the model file of the NMSSM as the 


input of the micrOMEGAs and set & = 0. 


e Since in the nMSSM, the LSP is too light to annihilate into Higgs pairs, there is no 


need to reconstruct the effective Higgs potential to calculate pem the annihilation 


channel y?X? > SS with S denoting any of Higgs bosons|53]. 


We thank the authors of the NMSSM Tools for helpful discussion on this issue when we finish 


such extension [19]. 


V. NUMERICAL RESULTS AND DISCUSSIONS 


With the above constraints, we perform four independent random scans over the param- 


eter space of the Type-II 2HDM, the L2HDM, the nMSSM and the NMSSM respectively. 


We vary the parameters in following ranges: 


1Xtanf «80, — 2/2 < sina € V2/2, m, < 30 GeV, As < 4m, 
5 GeV < mj, 4, < 500 GeV, 316 GeV < my+ < 500 GeV 


for the Type-II 2HDM, 


1 < tanß «80, — 2/2 < sino € V2/2, m, < 30 GeV, As € 4m, 
5 GeV < m,,4,, € 500 GeV, | 92GeV < my+ < 500 GeV 


for the L2HDM, 


0.1 X Ax0.7, 1 € tan B € 80, 100 GeV € m4 < 1 TeV, 
50 GeV < per, Mı < 500 GeV, —1 TeV < A, € 1 TeV, 0 € mg < 200GeV 


for the nMSSM, and 


0.1 € A,« € 0.7, 1 € tan f < 80, 100 GeV € my < 1 TeV, 
50 GeV < peg, Mı < 500 GeV, — 100 GeV < A, < 100 GeV 


for the NMSSM. 
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FIG. 4: The scattering plot of the surviving samples projected on m, — tan 8 plane. 


In performing the scans, we note that for the nMSSM and the NMSSM, some constraints 
also rely on the gaugino masses and the soft breaking parameters in the squark sector and 
the slepton sector. Since these parameters affect little on the properties of a, we fix them to 
reduce the number of free parameters in our scan. For the squark sector, we adopt the m7'^* 
scenario which assumes that the soft mass parameters for the third generation squarks are 
degenerate: Mg, = My, = Mp, = 800 GeV, and that the trilinear couplings of the third 
generation squarks are also degenerate, A, = A, with X, = A, — u cot 8 = —2Mg,. For the 
slepton sector, we assume all the soft-breaking masses and trilinear parameters to be 100 
GeV. This setting is necessary for the nMSSM since this model is difficult to explain the 
muon anomalous moment at 2c level for heavy sleptons|19]. Finally, we assume the grand 
unification relation 3M,/5a, = M»5/as = M3/a3 for the gaugino masses with a; being fine 
structure constants of the different gauge group. 

With large number of random points in the scans, we finally get about 3000, 5000, 800 and 
3000 samples for the Type-II 2HDM, the L2HDM, the nMSSM and the NMSSM respectively 
which survive the constraints and satisfy m, < 30GeV. Analyzing the properties of the a 
indicates that for most of the surviving points in the nMSSM and the NMSSM, its dominant 
component is the singlet field (numerically speaking, cos04 > 0.7) so that its couplings to 
the SM fermions are suppressed|15, |19]. Our analysis also indicates that the main decay 
products of a are 77 for the L2HDM[20], bb (dominant) and 77 (subdominant) for the Type- 
II 2HDM, the nMSSM and the NMSSM, and in some rare cases, neutralino pairs in the 
nMSSM|II9J. 
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FIG. 6: Same as Fig] but for Z > Fra. 


In Figl4] we project the surviving samples on the m, — tan plane. This figure shows 
that the allowed range of tan B is from 8 to 20 in the Type-II 2HDM, and from 37 to 80 
in the L2HDM. Just as we introduced before, the lower bounds of tan come from the 
fact that we require the models to explain the muon anomalous moment, while the upper 
bound is due to we have imposed the constraint from the LEP process ete~ — bbS — 4b, 
which have limited the upper reach of the bbS coupling for light S (for the dependence 
of bbS coupling on tan f, see Sec. II). This figure also indicates that for the nMSSM and 
the NMSSM, tan is upper bounded by 10. For the nMSSM, this is because large tan 8 
can suppress the dark matter mass to make its annihilation difficult (see and also Sec. 
II), but for the NMSSM, this is because we choose a light slepton mass so that large tan 3 
can enhance a, too significantly to be experimentally unacceptable. We checked that for 


the slepton mass as heavy as 300GeV, tan 6 > 25 is still allowed for the NMSSM. 
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FIG. 7: Same as Fig] but for Z — ay. 


In Fig] and Figl6] we show the branching ratios of Z — bba and Z — Tra respectively. 
Fig [5] indicates, among the four models, the Type-II 2HDM predicts the largest ratio for 
Z — bba with its value varying from 5 x 1079 to 6 x 1075. The underlying reason is in the 
Type-II 2HDM, the bba coupling is enhanced by tan 8 (see Fig), while in the other three 
model, the coupling is suppressed either by cot 6 or by the singlet component of the a. Figl6] 
shows that the L2HDM predicts the largest rate for Z + 7ra with its value reaching 1074 
in optimum case, and for the other three models, the ratio of Z — Tra is at least about 
one order smaller than that of Z — bba. This feature can be easily understood from the 
TTG coupling introduced in Sect. II. Here we emphasize that, if the nature prefers a light a, 
Z — bba and/or Z — Fra in the Type-II 2HDM and the L2HDM will be observable at the 
GigaZ. Then by the rates of the two decays, one can determine whether the Type-II 2HDM 
or the L2HDM is the right theory. On the other hand, if both decays are observed with 
small rates or fail to be observed, the singlet extensions of the MSSM are favored. 

In Fig[7] we show the rate of Z — ay as the function of ma. This figure indicates that 
the branching ratio of Z ay can reach 9 x 107°, 6 x 10719, 9 x 1071! and 4 x 10^? for the 
optimal cases of the Type-II 2HDM, the L2HDM, the nMSSM and the NMSSM respectively, 
which implies that the decay Z — ay will never be observable at the GigaZ if the studied 
model is chosen by nature. The reason for the smallness is, as we pointed out before, that 
the decay Z — ay proceeds only at loop level. 

Comparing the optimum cases of the Type-II 2HDM, the nMSSM and the NMSSM shown 
in Fig.5-7, one may find that the relation Brogpy > Brnussm > Brnmssm holds for any of 


the decays. This is because the decays are all induced by the Yukawa couplings with similar 
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FIG. 8: Same as Fig[5] but for Z > aaa. 


structure for the models. In the supersymmetric models, the large singlet component of 
the light a is to suppress the Yukawa couplings, and the a in the nMSSM has more singlet 
component than that in the NMSSM. 

Next we consider the decay Z — aaa, which, unlike the above decays, depends on 
the Higgs self interactions. In Fig] we plot its rate as a function of m, and this figure 
indicates that the Br(Z — aaa) may be the largest among the ratios of the exotic Z decays, 
reaching 107? in the optimum cases of the Type-II 2HDM, the L2HDM and the nMSSM. The 
underlying reason is, in some cases, the intermediate state h in FigB] (a) may be on-shell. 
In fact, we find this is one of the main differences between the nMSSM and the NMSSM, 
that is, in the nMSSM, h in Fig[3](a) may be on-shell (corresponds to the points with large 
Br(Z — aaa)) while in the NMSSM, this seems impossible. So we conclude that the decay 
Z — aaa may serve as an alternative channel to test new physics models, especially it may 
be used to distinguish the nMSSM from the NMSSM if the supersymmetry is found at the 
LHC and the Z — aaa is observed at the GigaZ with large rate. 

Before we end our discussion, we note that in the NMSSM, the Higgs boson a may be 
lighter than 1GeV without conflicting with low energy data from Y decays and the other 
observables (see Figl4l8). In this case, a is axion-like as pointed out in ld We checked 
that, among the rare Z decays discussed in this paper, the largest branching ratio comes 
from Z — abb which can reach 1.9 x 1079. Since in this case, the decay product of a is 
highly collinear muon pair, detecting the decay Z — abb may need some knowledge about 


detectors, which is beyond our discussion. 
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VI. CONCLUSION 


In this paper, we studied the rare Z-decays Z — f fa (f = b, T), Z > ay and Z — aaa in 
the Type-II 2HDM, lepton-specific 2HDM, nMSSM and NMSSM, which predict a light CP- 
odd Higgs boson a. In the parameter space allowed by current experiments, the branching 
ratio can be as large as 1074 for Z — f fa, 107? for Z — ay and 107? for Z — aaa, which 
implies that the decays Z > ffa and Z — aaa may be accessible at the GigaZ option. 
Since different models predict different size of branching ratios, these decays can be used to 


distinguish different model through the measurement of these rare decays. 
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